Determination of neutron-skin thickness using configurational information entropy* 
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Configurational information entropy (CIE) theory was employed to determine the neutron skin thickness of 
neutron-rich calcium isotopes. The nuclear density distributions and fragment cross-sections in 350 MeV/u 
40-60 Ca + °Be projectile fragmentation reactions were calculated using a modified statistical abrasion-ablation 
model. CIE quantities were determined from the nuclear density, isotopic, mass, and charge distributions. 
The linear correlations between the CIE determined using the isotopic, mass, and charge distributions and the 
neutron skin thickness of the projectile nucleus show that CIE provides new methods to extract the neutron skin 


thickness of neutron-rich nuclei. 
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I. INTRODUCTION 


Next generation radioactive nuclear beam facilities will 
provide new opportunities to explore extreme nuclei near and 
beyond drip lines. Nuclei with a large neutron excess can 
form exotic neutron skin or halo structures, which have at- 
tracted significant interest experimentally and theoretically 
for the past 30 years. The neutron skin thickness is de- 
fined as dnp = ôn — dp, which denotes the difference be- 
tween the point neutron and point proton root-mean-square 
(RMS) radii of a nucleus. Many methods have been devel- 
oped to experimentally determine the neutron skin thickness. 
However, most of the models are indirect measurements and 
model dependent. Typical methods used to determine neu- 
tron skin thickness include the reaction cross section (dR), 
charge-changing cross section (occ) [1, 2], electric dipole po- 
larizability [3], photon multiplicity [4], the m7 /m* ratio or 
x7 /=* [5, 6], 3H/%He ratio [7], and a decay half-life time 
[8]. The projectile fragmentation reaction, which is the main 
experimental approach for studying rare isotopes, is suitable 
for determining the neutron skin thickness owing to the ob- 
vious experimental phenomena induced by the neutron skin 
structure [9, 10]. For example, isospin effects in the isotopic 
cross section[11], neutron-abrasion cross section (Gpabr) [12], 
neutron removal cross section [13], mirror nuclei ratio or iso- 
baric ratio [14], and isoscaling parameter (œ) [15]. Parity- 
violating electron scattering (PVS) is the only method used to 
determine neutron skin thickness that is model-independent. 
Ref. [16] reports a theoretical investigate of PVS for 48Ca and 
208Pp, and a theoretical investigation of Bayesian approach is 
given in Ref. [17]. Determining the neutron skin thickness 
of 48Ca and ?°8Pb is presently of significant interest and is 
listed in the U.S. 2015 Long Range Plan for Nuclear Science 
[18]. The lead radius experiment (PREX) has been previ- 
ously used to determine the neutron skin thickness of 2°5Pb 
[19]. A recent PREX result indicates a much thicker neutron 
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skin than previously predicted [20]. Determining the neutron 
skin thickness of nuclei near the neutron drip line remains an 
important research outcome and is one of the most interesting 
topics in the new era of radioactive beam facilities. 


Information entropy theory was established by C.E. Shan- 
non [21]. The theory makes it possible to transform variables 
in a system into an exact information quantity [22] and has 
been used in various applications [23, 24]. The first applica- 
tion of information entropy theory in heavy-ion reactions can 
be traced to the study of nuclear liquid-gas transition in nu- 
clear multifragmentation [25]. Recent studies have extended 
it to study the information entropy carried by a single frag- 
ment produced in projectile fragmentation reactions, and re- 
vealed the scaling phenomenon of fragments covering a wide 
range of neutron excess [26-28]. Configurational informa- 
tion entropy (CIE) was developed to quantify the information 
entropy of a physical distribution [29], which connects the dy- 
namical and informational contents of a physical system with 
localized configurations. Many applications of CIE methods 
can be found in Korteweg-de Vries (KdV) solitions, compact 
astrophysical systems, and scalar glueballs (see a brief intro- 
duction in Ref. [23]), theoretical research of new Higgs bo- 
son decay channels [30], deploying heavier eta meson states 
in AdS/QCD [31], confinement/deconfinement transition in 
QCD [32], quarkonium in a finite density plasma [33], time 
evolution in physical systems [34, 35], etc. In projectile frag- 
mentation reactions, fragment distributions show a sensitive 
dependence on the change in neutron density [36-38], which 
makes it possible to determine the neutron skin thickness 
of neutron-rich nuclei. In this study, the CIE method was 
adopted to quantify the CIE of nuclear density and fragment 
distributions in projectile fragmentation reactions. The ana- 
lyzed data were generated using a modified statistical abra- 
sion ablation (SAA) model, which is known to be a good 
model for describing the fragment cross sections of projec- 
tile fragmentation reactions [39, 40]. 


Il. THEORIES 
A. Modified Statistical Abrasion-Ablation Model 


The modified statistical abrasion-ablation (SAA) model 
[39, 40] can be used for projectile fragmentation reactions 
at both intermediate and high energies, which improves the 
original SAA model by Brohm and Schmidt [41]. In quasi- 
free nucleon-nucleon collisions, the reaction is described as 
a two-step process: In the initial stage, the nucleons are de- 
scribed by a Glauber-type model as “participants” and “spec- 
tators”, where the participants interact strongly in an overlap- 
ping region between the projectile and target, while the spec- 
tators move virtually undisturbed [42]. In the second stage, 
the excitation energy is compared to the separation energies 
of protons, neutrons, and a to determine the type of particle 
the prefragment can emit according to min(Sp, Sn, Sq). After 
the de-excitation calculation, the cross sections of final frag- 
ments that are comparable to the measured fragments were 
obtained. The description of the modified SAA model is pre- 
sented in Refs. [23, 39, 40]. The colliding nuclei are com- 
posed of many parallel tubes oriented along the beam direc- 
tion. Their transverse motion is neglected, and the interac- 
tions between the tube pairs are independent. For a specific 
pair of interacting tubes, the absorption of the projectile neu- 
trons and protons is assumed to be in a binomial distribution. 
At a given impact parameter b, the transmission probabilities 
of neutrons (protons) of an infinitesimal tube in the projectile 
are calculated using 


ti(s — b) = exp{—[D,(s b)opil}, 


where DT is the normalized integrated nuclear density dis- 
tribution of the target along the beam direction for protons 
J sD} = Z" and neutrons f d?sD} = NT (N7 and Z7 are 
the neutron and proton numbers of the target, respectively). s 
and b are defined in a plane perpendicular to the beam direc- 
tion, and g; denotes the free-space nucleon-nucleon cross 
sections (2’,2 = n for neutrons and 7’, i = p for protons) [44]. 
The average absorbed mass in the infinitesimal tube limit at a 
given b is 


b)oni H Di (s 


(AA(b)) = f a*sDI(6)[1 — ta(s — b) 


+ J d?sDF(s)[1 — tp(s — b)]. (2) 


For a specific fragment, the production cross section can be 
calculated using 


o(AN, AZ) = I d?bP(AN, b)P(AZ,b), (3) 


where P(AN,b) and P(AZ, b) are the probability distribu- 
tions of the abraded neutrons and protons at a given impact 
parameter b, respectively. o(AN, AZ) is the residual frag- 
ment after the abrasion stage, which is called the prefrag- 
ment. The excitation energy of the prefragment is calculated 
as E* = 13.3(AA(b)) MeV, where (AA(b)) is the num- 
ber of abraded nucleons from the projectile, and 13.3 MeV 


is the mean excitation energy owing to an abraded nucleon 
[43]. In the second stage, the excitation energy is compared 
to the separation energies of protons, neutrons, and a to de- 
termine the type of particle the prefragment can emit accord- 
ing to Min(Sp, Sn, Sq). After the de-excitation calculation, the 
cross sections of final fragments that were comparable to the 
measured fragments were obtained. 

Fermi-type density distributions were adopted for protons 
and neutrons in a nucleus, as shown in the equation below: 


pe 
pi(r) = ——~__z, t=n,p (4) 
1+ exp( £74) 
where p? is the normalization constant of neutrons (2 = n) or 
protons (i = p), t; is the diffuseness parameter, and C; is half 
the density radius of the neutron or proton density distribu- 
tion. 


B. Configurational Information Entropy Method 


To determine the quantity of CIE incorporated in the frag- 
ment distributions, definitions of CIE were introduced. For 
a system with spatially localized clusters, when performing 
the CIE analysis, a set of functions f(x) € L?(R) and their 
Fourier transforms F (k) obey Plancherel’s theorem [45]: 


i] F(o)Pae = | |F()Pak, (5) 


where f(x) is square-integrable-bounded. The model fraction 
f(k) is defined as: 


|F (k)|? 
PIECE) Path’ 
where the integration is over all k, F'(k) is defined, and d is 
the number of spatial dimensions. 

The model fraction f(k) measures the relative weight of 


a given mode k. The quantity of CIE Sc|f] is defined as a 
summation of the Shannon information entropy of f(k) [21]: 


f(k) = (6) 


k 
= nia): (7) 
m=1 
Thus, the quantity of CIE has information about configura- 
tions compatible with certain constraints of a given physi- 
cal system. If all the modes k have the same mass, then 
fm = 1/N. The discrete configuration entropy reaches a 
maximum at Sc = ln N. If there is only one mode, Sc = 0 
Continuous CIE can also be defined for continuous distri- 
butions, such as the nuclear density distribution. For non- 
periodic functions in the interval (a, b), 


- | FE mif aak, 8) 
where f(k) = f(k)/f(k)max [ J (k)max is the maximum frac- 


tion]. The normalized function f(k) guarantees that f(k) < 
1 for all k modes, and f(k) ln f(k) denotes the CIE density. 
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Fig. 1. (Color online) Fermi-type nuclear density distributions [in 
panels (aj)] and the corresponding FFT spectrum [in panels (b7)]. 
Panels (al), (a2), and (a3) show the proton densities (pp), neutron 
densities (pn), and the nuclear density difference Ap = pn — pp for 
40-60 Ca isotopes. 


IHI. RESULTS AND DISCUSSION 


The 350 MeV/u 4?Ca + °Be reactions were calculated us- 
ing the modified SAA model (A, refers to even mass numbers 
from 40 to 60). The cross sections of fragments with Z rang- 
ing from 3 to 20 were obtained. For the sake of clarity, only 
part of the calculated results are shown in the figures. Fig- 
ure | is a plot of the Fermi-type nuclear density distributions 
and their fast Fourier transformation (FFT) spectra. An ob- 
vious increase in pn is observed from 4°Ca to ®°Ca, whereas 
the opposite trend is observed for pp. A two-peak structure 
is evident in the FFT spectra, where the second peak is lower 
than the first. The difference between the neutron and proton 
density distributions Ap = pn — pp is also shown. For 40Ca, 
Ap is very small, while Ap increases as the neutrons in the 
projectile increase. The peaks in the FFT spectra of p, and 
Pp are not clearly shown. Based on the FFT spectra f(k), the 
CIE of pn, Pp, and Ap can be determined using Eq. (7), which 
are denoted by S£ [f], SE [f], and SÊ” [f], respectively. 

The isotopic cross section (az) distributions produced in 
the 350 MeV/u 4?Ca + Be reactions are plotted in Fig. 2. In 
panels (ai), from Zf = 7 to 20, the isotopic cross-section dis- 
tributions in the 4°~®°Ca reactions are similar for fragments 
with small Z;,, while a shift to the neutron-rich side is ob- 
served for larger Z. The symmetric Guassian-like shape 
of the isotopic distribution is altered by the enhanced cross 
sections of neutron-rich fragments in neutron-rich reaction 
systems, showing the isospin effect in fragment production 
induced by the increased neutron density on the surface of 
neutron-rich nuclei [40]. The FFT spectra of the isotopic dis- 
tributions are shown in Fig. 2 (bj). In each FFT spectrum, 
only one peak is observed. The amplitudes of the FFT spec- 
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Fig. 2. (Color online) The isotopic cross section distributions of 
fragments in the 350 MeV/u 4>Ca + °Be reactions calculated using 
the modified SAA model. Panels (a7) correspond to Z = 7, 12, 16, 
and 20 isotopes, respectively, and panels (bj) plot the corresponding 
FFT spectra of (ai). 


tra of different isotopic distributions decrease as the projectile 
becomes more neutron-rich, except for Z = 20. Based on the 
FFT spectra of oz distributions, the quantities of CIE of the 
isotopic distributions are determined according to Eq. (7), 
which is denoted by S27 [f]. 

The correlation between the Sc[f] of the density distribu- 
tion and np of the projectile nucleus is shown in Fig. 3(a). 
Both S£ [f] and S [f] decrease linearly with an increase in 


Sap from 4°Ca to ©°Ca. S&?[f] also decreases linearly with 
increasing Opp, except for *°Ca. The correlation between the 
Sé7(f] of different Z and dnp of the projectile nuclei are 
plotted in panel (b). The S27 [f] of isotopes from Ze = 10 
to 18 are also found to decrease with the increasing dnp of 
the projectile nucleus. The 527[f] of the fragment near the 
projectile nucleus was more sensitive to the change in dnp. 


The mass yield (ø 4) distributions in the 350A MeV ArCa + 
°Be reactions are shown in Fig. 4. In each reaction, the mass 
yield increases with the As, of the fragment until it is close 
to the projectile nucleus. In different reactions, a very simi- 
lar trend of mass distribution was observed, which decreased 
with the increasing mass number of the projectile nucleus. 
The corresponding quantities of CIE were determined from 
o distributions, which are labeled as S2“[f]. The correla- 
tion between S¢^ [f] and ônp for projectile nuclei is shown 
in Fig. 4 (b). Except for the bend point formed at dnp for 
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Fig. 3. (Color online) (a) The correlation between the CIE of pn 
(squares), pp (circles), Ap (triangles) and neutron skin-thickness dnp 
of “°Ca. (b) Correlation between the CIE of the isotopic cross- 
section distributions S27 [f] and dup of “PCa. The lines denote the 
linear fitting of the correlations. 


42Ca owing to the transition of proton-skin to neutron skin, 
the S2[f] ~ dnp correlation was found to be linear for reac- 
tions of Ap > 44. 

The charge cross section is defined as the summation of the 
isotopic cross sections øc = )>,, ao (Ao, Z). The charge cross 
section distributions in the 3504 MeV “Ca + °Be reactions 
are shown in Fig. 5. Similar trends of øc distribution trends 
similar to those of g a were observed. The determined CIE of 
oc distributions, labeled as Sé°[f], were linearly correlated 
to the neutron skin thickness of the projectile nuclei. 

The CIE approach transfers the experimental distributions 
to quantified parameters and provides information probes for 
determining the properties of a system. From the S27[f] ~ 
Onp, Sc lf] ~ Onp, and Sé°[f] ~ dnp correlations, it was 
observed that the CIE determined from isotopic, mass, and 
charge distributions decrease with increasing neutron skin 
thickness, respectively, and they had good linear correlations. 
The determination of neutron skin thickness, in particular, nu- 
clei near the neutron drip line, is limited by the unavailability 
of effective probes. The linear correlation between the CIE 
and neutron-skin thickness of neutron-rich nuclei provides 
new approaches to determine the neutron skin thickness of the 
projectile nucleus by measuring the fragment distributions in 
projectile fragmentation reactions. 
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Fig. 4. (Color online) Panel (a): The mass yield distributions of 
the 350 MeV/u 4?Ca + °Be reactions calculated using the modified 
SAA model. Panel (b): The correlation between S24 [f] determined 
from the mass yield distribution and the np of “PCa (Ap denotes 
even mass numbers from 40 to 60). 


IV. SUMMARY 


With the vast opportunities for very asymmetric nuclei 
available in the new era of radioactive ion beam facili- 
ties, neutron skin thickness is one of the most important 
questions in nuclear physics. In this study, CIE theory is 
adopted to quantify the information entropy incorporated in 
nuclear density distributions and fragment cross-section dis- 
tributions in 350 MeV/u 4°-©°Ca + °Be projectile fragmen- 
tation reactions calculated using the modified SAA model. 
CIE quantities of nuclear density distributions (S£ [f] and 
see [f]), isotopic cross-section distributions (S27 [f]), mass 
cross-section distributions (S¢^ [f]), and charge cross-section 
distributions (S2°|f]) were determined. The correlations be- 
tween SRS] ~ Sup, SELF] ~ Sup, SEPI] ~ Snp, and 
SE [Sf] ~ Onp, SC lS] ~ dnp, and SEC [Ff] ~ dnp were 
also investigated. For neutron-rich calcium projectiles, ob- 
vious linear dependences of S® [f], S> [f], and SÊ” [f] on 
dnp Were observed. The Sg*7 [f] of fragments with different 
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(Color online) Similar to Fig. 4 but for the charge cross 


section distributions. 
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